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Abstract. A conserved lattice gas witli random neighbor hopping of active 
particles is introduced which exhibits a continuous phase transition from an active 
state to an absorbing non-active state. Since the randomness of the particle 
hopping breaks long range spatial correlations our model mimics the mean-field 
scaling behavior of the recently introduced new universality class of absorbing 
phase transitions with a conserved field. The critical exponent of the order 
parameter is derived within a simple approximation. The results are compared 
with those of simulations and field theoretical approaches. 



Recently Rossi et al introduced a conserved lattice gas (CLG) with a stochastic 
short range interaction that exhibits a continuous phase transition to an absorbing 
state at a critical value of the particle density |^]. The CLG model is expected to 
belong to the universality class of absorbing phase transitions with a conserved field . 
This universality class is different from the well known universality class of directed 
percolation (see ^, ^, or for an overview Q). In this work we introduce a modified 
CLG model with a random neighbor interaction. This random neighbor interaction 
suppresses long range correlations and the model is expected to be characterized by the 
mean field scaling behavior. The critical exponent of the absorbing phase transition is 
derived within a simple approximation. Numerical simulations of various lattice types 
with different numbers of next neighbors confirm the obtained results. 

Let us consider a system consisting of L sites on a chain with periodic boundary 
conditions. At the beginning one distributes randomly TV = pL particles on the 
system where p denotes the particle density. A particle is called active if at least one 
of its two neighboring sites is occupied. In the original CLG model active particles 
jump in the next update step to one of their empty nearest neighbor site, selected at 
random 0. In the steady state the system is characterized by the density of active 
sites Pa which depends on p. The density of inactive sites is given by p — pa and 1 — p 
is the density of empty sites. The density pa is the order parameter of the absorbing 
phase transition i.e., it vanishes at the so-called critical density of particles pc 
(which of course is lower than the trivial value p = 1/2). 

In our modification of the CLG model active particles are moved to a randomly 
chosen empty lattice site which suppresses long range correlations. On the other hand 
short range correlations still do exist. For instance close to the critical point, where 
the density of active sites is sufficiently low, the density-density correlation function 
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displays an alternating behavior due to the repulsive interaction of the particles. But 
these correlations are short ranged and it is therefore possible to neglect them. The 
values of the critical exponents are not affected by this approximation. 

Let us consider a configuration C of the lattice with n active particles. The 
number of active sites may change in each particle hopping, i.e., in each elementar 
update step. For instance if both new neighbors of the hopped particle are empty 
the number of active particles is reduced by one. An = —1. Without correlations 
the corresponding probability is (1 — p)^. If one of the new neighbors of the hopped 
particle is occupied by an inactive particle (j> — p ~ Pa) and the second neighbor is 
empty (p = 1 — p), the number of active sites is increased by one (An = 1), and 
the corresponding probability is given by p = 2(1 — p){p — pa)- AH other possible 
configurations and the corresponding probabilities are listed in table |l|. 

In this way one can calculate the probabilities that the number of active particles 
are changed by An and one gets 

Pa„=-i - (1-p)' 

Pa„=o = 2pa(l -p) + p^ 

(1) 

Pa„=i = 2 (p - Pa) (1 - p) + 2pa (p - Pa) 
Pa„=2 = (P - Pa)^ 

The expectation value of An is 

2 

£;[An] = AnpA„ ^ -l-2pa + 4p-p^ (2) 

An=-1 

The average number of active sites is constant in the stationary state, i.e., the 
expectation value of An should be zero in the steady state. Using the constraint 
£'[An] = it is possible to calculate pa as a function of p and one gets 

4p^p'-l 

Pa = ^ ■ (3) 

The corresponding curve is plotted in figure|l| for < p < 1. Negative values of 
Pa corresponds to an absorbing state (i.e. pa = 0). The critical point is determined 



Table 1. The configuration of the lattice before (C) and after (C') a particle 
hopping. Only the target lattice site where a particle hops and its left and right 
neighboring sites are shown. Empty sites are marked by o, inactive sites are 
marked by *, and active sites by •. An denotes the change of the number of 
active sites due to the particle hopping and p is the corresponding probability of 
the configuration C if one neglects spatial correlations. 
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Figure 1. Figure The density of active sites pa as a function of the particle 
density p in the steady state. Negative values (dash line) correspond to an 
absorbing state with pa = 0. 



by Pa = 0, which leads to pc 2 — \/3 (the second sohition can be neglected since 
pc = 2 + \/3 > 1). Writing pa as a function of the reduced density Sp = p — pc one 
gets 



(4) 



Thus the order parameter of the absorbing phase transition vanishes in leading order 
as Pa Sp, i.e., the order parameter exponent is /3 = 1. 

We briefly remark that it is straight forward to generalize the above derivation 
from a chain with two neighbors to a c?-dimensional cubic lattice with z — 2d 
neighbors. Since An depends on the number of inactive sites (*) only (see table|l|) the 
corresponding probabilities are just polynomials in = p — pa, i.e., 

pa„=-i = (l-p)^ 
pa„=o = (1-p,)- - (l-p)^ 

Pa„>. = (Ajpf"(l-P*)^"^"- 

Using again the steady state condition E[An] = one gets 

Pe. = p ■ 



(5) 



(6) 

The critical density pc is determined by pa ~ 0. Expanding equation (^) around pc 
yields 

Pa = (2-pc)5p + 0(V), (7) 

i.e., the order parameter exponent is /3 = 1 independent of the number of next 
neighbors. 

In the following we compare this results with those obtained from simulations. 
We simulated a one dimensional chain (z = 2) and a two dimensional square lattice 
[z = 4) where the active particles are moved to a randomly chosen empty lattice site. 
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Figure 2. Figure The average density of active sites for a one-dimensional cliain 
of size L with random particle hopping = 2). The dashed lines corresponds 
to a linear fit, i.e., the order parameter exponent is /J = 1. The inset displays 
the fluctuations of the order parameter which exhibits a discontinuous behavior 
at the critical density pc (dash line). 



In both cases a random sequential update was used and the results are plotted in 
figure^ and figure ||, respectively. As expected the critical density decreases with 
the number of next neighbors z. Note that the critical value of the chain « QflTl^ 
and of the square lattice pc ~ 0.1244 differs slightly from the above analytical result 
Pc — 0.2679... and pc = 0.1380..., respectively. This deviation, which decreases with 
increasing z, is caused by the neglection of correlations between neighboring sites. 
Simulations reveal that these correlations exists but are of short range (not shown). 
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Figure 3. Figure The density of active sites for a two-dimensional square lattice of 
linear size L with random particle hopping (2: = 4). The dashed lines corresponds 
to a linear fit, i.e., the order parameter exponent is /3 = 1. The inset displays 
the fluctuations of the order parameter which exhibits a discontinuous behavior 
at the critical density pc (dash line). 
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Therefore, the critical value is shifted but the scaling behavior itself agrees with our 
analytical results, i.e., the order parameter vanishes at the critical point continuously 
with an exponent /3 = 1. 

Additionally to the order parameter its fluctuations 



are measured in the simulations. Here D denotes the dimension of the system. As can 
be seen from the insets of figure^ and figure^ the fluctuations exhibit a discontinuous 
behavior (jump) at p^. 

The value (3=1 and the jump of the fluctuations was also observed in the CLG 
model above the critical dimension Z?c = 4 |^ where the scaling behavior of the model 
is determined by the mean-field exponents. Furthermore the values /3 = 1 and = 4 
were predicted within a field theoretical approach which is expected to represent the 
universality class of absorbing phase transitions with a conserved field . 

Finally we just mention that the derivation of our results corresponds to a 
mapping of the dynamics to a branching process (see for instance [||). There each 
active site can create i S {0,1,2,3} active sites in the next generation with the 
probability pi = PAn=i-i- The steady state condition corresponds to the condition 
that the average number of created active sites in the next generation is one which 
yields directly equation (^) . 



S. L. would like to thank H. K. Janssen for helpful discussions and useful comments 
on the manuscript. 
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